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Abstract 



We calculate exactly the modes of motion of the Time-averaged Orbiting 
Potential (TOP) trap with its four degrees of freedom, namely the three 
translations and the spin, taken into account. We find that, when gravity 
is neglected, there are two parameters in the problem namely, the angular 
velocity of the rotating field and its strength. We present the stability di- 
agram in these parameters. We find the mode frequencies calculated from 
the time-averaged potential model used by the inventors of the TOP is an 
excellent approximation to our exact results. However, for other parameters, 
this may not be the case. 



1 Introduction. 



The first observation of Bose-Einstein condensate (BEC) ||T| was done in 
tlie so-called Time-averaged Orbiting Potential (TOP) trap. Although later 
other magnetic traps, optical plug, and even an all-optical one0] were used, 
the ingenious TOP trap continues to be a workhorse in the trade P] because 
of its unique advantage as pointed out by |Q. It seemed to us thus worthwhile 
to extend the analysis given in Ref.0], which is based on the concept of a 
time-averaged adiabatic potential, by studying the exact motion of all the 
degrees of freedom in the problem, namely the three translational degrees 
of freedom and the spin degree of freedom of the particle 0. This has the 
advantage of allowing the calculation of limitations required by stability. 
We find that, when gravity g is neglected, there are two parameters in the 
problem, the angular velocity of the rotating field Q and its strength a, and 
we present the stability diagram in the a-fl plane. The parameters reported 
in Ref.[0 are shown to be well inside the stable region. Their calculated 
frequencies agree excellently with our exact frequencies for their experiment. 
When the limit of high field strength is taken, our analytic results reduce to 
their formulae. The stability diagram that we found shows that the TOP 
trap is very flexible for the experimentalist in terms of allowed parameters a 
and fl. Our treatment is classical, but we also discuss what is involved in a 
quantum-mechanical calculation. 

The structure of this paper is as follows: In Sec.0 we first define the 
problem, write down the equations of motion and find stationary solutions 
of these equations. Next, we perturb the stationary solutions, linearize the 
equations of motion and derive the secular equation which determines the 
frequencies of the various possible modes. The secular equation is given as a 
function of a, Q and g. In Sec^ we use the secular equation to find the mode 
frequencies for the TOP trap described in Ref. |Q, and compare our exact 
results with their calculations and measurements. In Sec.| we specify to the 
case g = and derive approximate expressions for the mode frequencies in 
the limit where a is large, starting from the secular equation. These are 
found to agree with the corresponding expressions given in Ref. [Q. In Sec.|^ 
we present the stability region in the a-Q plane (for g = 0) for which stable 
trapping occurs, and comment on its implications and limitations. Finally, 
in Sec.|^ we summarize our results and discuss briefly the related quantum- 
mechanical problem. 
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2 Derivation of secular equation. 



2.1 Mathematical formulation and physical interpre- 
tation. 

We consider a particle of mass m, intrinsic spin S and magnetic moment 
/i = —fih, moving in 3D space in the presence of uniform gravitation —Gz 
and an inhomogeneous time-dependent magnetic field given by 

H — H' (^~2PP + (cos (v? — Qrt) P — sin {(p — flri) 'p) ■ (1) 

Here, (/i, p, </?) are the height, radial distance and polar angle of the center 
of mass of the particle in cylindrical coordinate system, (z, p, (p) are the cor- 
responding unit vectors, H is the strength of the (uniform) time- dependent 
field rotating at an angular velocity Vt^ about the ^-axis, and H' is the time- 
independent field gradient in the 2;-direction. 

The equations of motion for the center of mass of the particle r = zz+pp 

are 

= -pV (n • H) - mGz, (2) 
and the evolution of its spin is determined by 

= -ph X H. (3) 

(Jjv 

In the following, it is convenient to express n- a unit vector in the direction of 
the spin, in terms of its components along the (z, p, ip) directions. We denote 
these by n^, Up and n^, respectively. Note however, that the unit vectors 
(p, (f) themselves depend on time according to 

d(p _ _d(p^ 
dt ~ dt^' 
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Substitution of Eq.(|T]) into Eqs.(|^) and (j^), and making use of Eq.(^ yields 



fiH' 
2m 



-Ur. 



dp d(p d (fi pH' 
d^h pH' 



dn 



dn 



_p _d^ 
dt dt""^ 



dt^ 



G 



m 



_^ dip 
dt ^ dt""' 



p 

's 



drty 



— — [H'hn^p + Hriz sin — Vtrt)] 
Hriz cos {if — Qrt) — -pUzH' — UpH'h 



p 



[—n^H cos {if — fifi) + ^pn^H' — UpH sin {ip — fi^t)] 



(5) 



dt S 

Eqs.(l^) are almost the equations of motion of the particle in a coordinate 
system which is rotating with the field. The only difference is in the definition 
of the angle <y9, which is measured with respect to the fixed x axis rather than 
with respect to the axis defined by the rotating field. To show this we rewrite 
the equations of motion in the rotating frame by substituting 



dt^ 



dr dr ^ ^ 

, Tt Tt^^'^ " ^' 

— - + 2f2rZ X h f2^zx (z X r) 

it^ , , dt 



du\ 



dt I 



dh 



— ] ^ [ — ] + ^IrZ X n. 



dt 



This brings the equations of motion Eq.(|^) and @ into 



m 



and 



d'^v dv 

+ 2f2rZ X — + f2^zx (z X r) 



dt^ 



dt 



dt 



-pnx 



-pV 



n H 



S 



H f^^z 



(6) 



mGz 



where now the magnetic field H is time-independent and is given by 



H =H' 



-pp + hz ) +ifx. 



(7) 
(8) 

(9) 
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Eqs.(l^) and (j^) indicate that in the comoving frame, the particle is acted 
upon by an additional uniform magnetic field — {Sflr/lJ') z, a centrifugal force 
—mQf.zx (z X r) and a velocity-dependent Coriolis force mv x (2fij.z). The 
action of the Coriolis force may also be interpreted as the Lorentz force of the 
inertial field ~ fi^z which acts on the mass of the particle. Finally, rewriting 
Eqs.(l^) and (|]) in polar coordinates yields Eqs.(|) with ip — Qrt replaced by 

As the number of parameters in the problem is relatively large, we rewrite 
the equations of motion in terms of normalized coordinates. We thus define 

/ ^2 \ 1/3 

Ro = 



\fimH' 

as the characteristic length-scale in the problem, and 

1/3 

/ (iiH'r \ 

^0 = 



mS 



as the characteristic angular velocity. This allows to define the dimensionless 
quantities 

r = p/Ro 
z = h/Ro 

T = Q^t 

g ^ G/{nlRo) 



with which Eqs.(|^) become 



(Pr 



dr J 

dr dip d (f 

2 h r 

dr dr dr'^ 



1 

2 ^ 
1 

-n,n 



dup dip 



dr^ 



(10) 



dn,, 



, i-n^, = —zriu, — ariz sin iip — Vtr) 

dr . dr 



dip 



driz 



dr ^ dr ^ 



—ariz cos {ip — Vtr) + hrriz + UpZ 



n,pa cos {ip — Qr) — -rn^, + Upa sin {ip — Qr) 
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In this form, one is left with only three parameters 

• a-the normalized strength of the rotating field, 

• fi-the normalized angular speed of the rotating field, and 

• g-the normalized free-fall acceleration. 



2.2 The stationary solutions and their stabiUty. 

We seek a solution in which the particle moves synchronously with the field 
at a constant radius and height. Setting 

(f = ipo + Qr ; r(r) = vq ; z{t) = zq 

in Eqs.([l0D we find two possible solutions, given by 



Up = — 2n^ro 
= 

n. = -g (11) 
ifo = 90° ±90° 

The value of vq is determined by the condition that |n| = ^ra^ + n"^ + = 1, 
giving 



1-r 

(12) 



■ 

It can be easily shown that the stationary solution corresponding to <^o = 
180° has its magnetic moment antiparallel to the direction of the local mag- 
netic field, whereas for the ipQ = 0° solution it is parallel to the direction of 
the field. From Eq.(|T^) we also conclude that for a stationary solution to 
exist, \g\ must be smaller than 1. This is simply a consequence of the fact 
that, in our model, the magnetic field cannot apply a force greater than fiH' 
in the +z-direction. When the weight of the particle mG is greater than 
fj^H', the magnetic force cannot balance the weight of the particle. The lat- 
ter then accelerates in the axial direction, and no stationary solution exists. 
Note also that the particle is located above the origin with its spin pointing 
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inward, even in the absence of gravity 0. This is necessary in order to have 
a 2;-component of the field. The latter exerts torque on the spin that causes 
it to rotate synchronously with the field. 

To check the stability of the solutions found, we add first-order pertur- 
bations. We set 



ZiT 



r(r) = To + 6r 
^ = fir + 90° ± 90° + 5ip 

Up = — 2f2^ro + Sup 
n^ = + 5n^ 
nz = -g + Sn^, 



+ 6z 



(13) 



substitute these into Eqs.(0), and retain only first-order terms. We find 
that the resulting equations for the perturbations are 



2firo 



dr 
d^z 



4fiVo 



6n, 



dSrip 
~d^ 



drj 
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(14) 



2fiVo Vr°^")^"-^"^^^ 



dSn, 



1 o (^'"O ± Oi) 



dr dr 
In deriving these equations we used the constraint 

h-6h =2Vl?rQ5np + gduz = 0, 



(15) 



since n is, by definition, a unit vector. We have also discarded in Eqs.(|l^ 
the equation corresponding to the last equation in Eqs.(0) as it turned out 
to be identical to the fourth. 

Looking back at Eqs. (|T^) we note that the two possible solutions differ by 
the sign of a (reversal of the direction of the rotating field). It is therefore 
suffices to concentrate on the solution with ipo = 180° for example, and to 
study both positive and negative values of a. This is what we do in the 
following. 
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To look for oscillatory (stable) solutions for Eqs. ( |T^ ) we set 



6r = ((5r)o e"^^* 
6z 



{6z) 
(5n, 



-iu)t 



-iujt 



2/0 



-ILdt 



<fiJO ' 



inside Eqs.(0) and get 
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(16) 



This equation has non-trivial solutions whenever the determinant of the ma- 
trix A vanishes. Thus, the equation 



where 



16rofi^ det A =A4UJ^ + Asto^ + Aaw^ + Aiu^ + Aq = 0, 



Ao = -32Q^^arl 



(17) 



(18) 



+ 2n^aro - 2g%n^a 



A2 = -lenv^ - 8a^n 



0^ - 8a^n^ 

2 Qon6„3^ in„n2„ 



A 



- 16fiV^ - 32nV^a - lOafiVo 

- 16n^ grl - 2^Vo^]^a - g'^r'^n^ 
3 = + + 4aro + 32^^^ 

4„2 



Aa = -i6n\ 



0' 



determines the eigenfrequencies u of the various possible modes. 
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At this point we pause and prove that the solution with (fQ = 0^ is not 
stable for any combination of a, Q and g: We consider the fourth order 
polynomial Eq.(|l^) as a polynomial in x = cu^. We note that when x = 0, 
the polynomial takes on the value Ao = — 32f2^°arQ which is positive when 
a is negative. When x — oo on the other hand, it takes on a value which 
is asymptotic to A^x^ = —IGQ^^r^x^ which, for sufficiently large (negative) 
X, gives a negative value. Thus, at least one root x = tu^ of the polynomial 
must be negative, corresponding to two purely imaginary frequencies with 
opposite signs, and hence to an unstable solution. 

3 Analysis of the TOP trap of Ref. [@|| . 

We now put Eq.(|T^) to use by calculating the mode frequencies of the TOP 
trap described in Ref. 0]. The parameters in this trap were: H = lO^^Tesla, 
G ~ 10m/ sec^, H' = 2 AT esl a /meter, m ~ {Aju,/ AH)mproton = 1.416- 
10^^^[ii'(?] (where Aju, is the atomic mass of Rubidium, Ah is the atomic 
mass of Hydrogen and nip is the mass of proton), f rotation = 7.5KHz, /i = 
= 4.6 ■ 10~'^'^ Joule /Tesla and S = h = 1 ■ 10~^'^ Joule ■ sec. From these 
parameters we find that flo = 2.049 ■ lO'^rad/sec, Rq = 1.856 ■ lO^^meter, 
a = 2.245 ■ 10^, g = 0.1283, = 2.3, ro = 0.0937, zq = -288, Up = -0.99, 
riz = —0.1283 and = 0. These correspond to a particle orbiting at a radius 
of po = 17.4?T,m and height zq = 53.5fim above the origin. The magnetic field 
at that point is about — (p— 0.128z) • lO^^Tesla, and is therefore dominated 
by the strength of the rotating field. 

Solving for the four roots u of the characteristic polynomial Eq. (p!7| ) and 
multiplying the tu's by Qq/2ti gives the frequencies: 7.38MHz, 7 .h2AAAlh22K H z , 
7.475558498iri7z and Q7mHz. The highest frequency 7.38MHz corre- 
spond to the precessional mode of motion. Its frequency is very close to 
the value ^Hs/2tiS where Hg is the strength of the local magnetic field at 
the stationary point. As was calculated above, Hg is roughly equal to the 
strength of the rotating field H. The two middle frequencies ~ 7.5KHz 
turn out to be very close to the frequency of rotation of the field 7.5KHz. 
The reason that we keep so many digits in these frequencies will be clar- 
ified shortly. Calculation of the eigenvectors for these two modes yields 
{Sr)^ / {Sz)q ~ 1.9 ■ 10^ ^ 1 for both, which shows that their coupling to 
the axial coordinate is small, indicating that these are the two lateral vi- 
brational modes. Recall that the frequencies that are found here correspond 
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to the rotating frame. To find the frequencies in the laboratory frame we 
exphcitly calculate the time dependence of the lateral coordinates: The x 
coordinate for example, is given by x (r) = {tq + 6r) cos [180° + Qt + 6ip] 
where 6r = {Sr)^ cos (a;r), 6ip = {Sip)^ cos {lot + 0o) and 0o is a phase which 
is determined by the eigenvectors. Expansion of x (r) to first order in the 
perturbations shows that it contains a term proportional to the product 
cos [180° + Qt + 6lp] cos {ujt + 0o) which contains a term with slow frequency 
u — Q and another term with fast frequency u + Q with the same amplitude. 
We thus find that, in the laboratory frame, the frequencies of the axial vi- 
brational modes are given by subtracting and adding the rotation frequency 
from the middle frequencies found by using Eq.(17).This gives the two slow 
lateral frequencies +2AAAlh22Hz and -24.441502//^, and the fast lateral 
frequencies 15.0244KHz and 14.9756KHz. The eigenvector for the lowest 
frequency 67.99Hz satisfy {Sr)^ / {Sz)q ~ 5 ■ 10~^ -C 1, which is a clear indi- 
cation that this is the axial vibrational mode. In this case, the frequency of 
vibration is the same both for the rotating frame and the laboratory frame 
so no subtraction/addition of Q is needed. 

Table I compares our results to the measured and calculated results re- 
ported in 1^. We conclude that the Time-averaged Orbiting Potential (TOP) 
approximation is indeed a very good approximation to the exact result ob- 
tained here. Note however, that the two fast lateral frequencies u + Q that 
our analysis yield, were not reported in Ref. possibly because they were 
not looked for. 

Note also that in the time-averaged orbiting potential approximation, 
the resulting potential is necessarily isotropic in the lateral plane. Hence, 
in this approximation there is only one lateral vibrational frequency. The 
exact analysis presented here, on the other hand, gives naturally two fre- 
quencies, corresponding to the two slow lateral vibrational modes, and two 
more frequencies, corresponding to the fast lateral modes. For the TOP trap 
parameters given above we find the slow frequencies to be very nearly equal, 
up to the seventh significant digit. The splitting is a consequence of the fact 
that the clock-wise lateral vibrational mode is not equivalent to the counter 
clock-wise lateral vibrational mode since in both cases the spin precesses 
in the same direction. This splitting may be large for other choice for the 
parameters. 
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4 Approximate expressions for the mode fre- 
quencies. 



For the numerical example of the previous section we have found that (7 <C 1, 
a S> 1 and ~ 1. These relations hold also for more recent TOP traps 
0, and it is therefore natural to find approximate expressions for the 
mode frequencies under these limits. We expect that for these values, the 
four frequencies will be as follows: A very high frequency ujp corresponding to 
the precession, two very close middle frequencies u^y which are almost equal 
to Vt and correspond to the lateral vibrations, and a very low frequency uj^ 
corresponding to the axial vibrations. The large dynamic range spanned 
by these frequencies suggests that there is a certain relation between the 
coefficients of the polynomial in the secular equation Eq. (p^Tf) . We exploit this 
relation in order to find approximate expressions for the mode frequencies: 
We first construct a fourth-order polynomial in w^, whose roots 0;^, uly and 
ojI satisfy the relations oj"^ ^ oj^y ^ u^, expand it in powers of o;^, and keep 
only dominant terms in the coefficients. This gives 



{u^ - ul) {u^ - ul/ {u^ - ul) (19) 

Comparing this with Eq.(|l^) shows that 

ul ^ hm (~^] = a' (20) 

oj^ ~ hm — — = — . 

5,=0,a>l \ AiJ a 

Clearly, this approximation is not sufficient to determine the difference of the 
frequencies of the lateral vibrational modes. To furnish these differences, we 
substitute a;^ — >• + d into Eq.(|17D and expand to second order in d. The 
result is 

d' X: ^^^^^n {^T' + dY,nA^ {^fy-' + ^ A„ = 0. (21) 

n=0 n=0 n=0 
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Setting g = and solving for d gives 



d = ±^ + 0{a-'), (22) 
V 2a 



corresponding to the slow laboratory-frame frequencies 

1 

In particular we find that 



);^, = V^^^Td -Q = ±^ + O {a-^/^) . (23) 



± V8 + O , (24) 



which shows that our analysis reduces to the results reported in Ref. 0]. 

In non-normalized units the mode frequencies in the laboratory-frame are 
given by 



1 1 1 ~ 



5 The stability region for ^ = 0. 

Having analyzed the TOP trap for the limit a ^ 1, we turn our attention to 
study the region a, ~ 1, which has not yet been exploited experimentally. 
To keep matters simple, we specify to the gravitation-free case g = 0. We 
use Eq.([T7|) to scan the a-fl plane in the search of points corresponding to 
stable solutions. The result is shown in Fig.(H). 

This figure shows the boundary lines between points corresponding to sta- 
ble regions and unstable regions. The fact that these lines consist of different 
segments, indicate that two different segments of a given line correspond to 
different modes that become unstable. For example, going upward along 
the a = 2.5 line, we find that the two slowest modes coalesce and become 
unstable at ^ 0.73. As Q is increased, these modes become stable again 
at ~ 0.91. When Q is further increased, the two fastest modes become 
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unstable at Q ~ 1.72. Note that we have found a small stability region in 
the range Q ^ 1.2-2 and a ^ 0-0.05, which we did not investigate in detail. 

As Aq oc Q^^a we conclude that points, both along the ^2 = line and 
along the a = line, have one mode with vanishing frequency, corresponding 
to a soft mode. In addition, all points in the a < half-plane correspond 
to unstable solutions, as was proved earlier. Note also that the coefficients 
in the secular equation really depend on and not just on Q. Hence, the 
continuation of the stability diagram to the 1] < half-plane is simply a 
mirror reflection of the f2 > half-plane with respect to the f2 = line. 

The stabihty diagram shows that the TOP trap is much more tolerant 
that what one would have expected. There are many points near the cor- 
ner of the first quadrant of the a-Q plane that may be used experimentally. 
Note however that as tq ~ 1/^^^, the use of too low a value for Q results 
in a large radius. This radius may fall outside the region in which our lin- 
ear approximation to the spatial dependence of the field holds. The lower 
bound on Q is therefore determined by the second derivative of the field. 
In addition, quantum-mechanical considerations set a lower bound on a as 
well. It can be shown that the extent of the wavefunction of the particle 
Axquantum ~ sjhjmu)^ and the extent of the field Axfieid ~ H/H' scale as 
Aa;quantum/Aa;fieid ^ oT"^!^. Thus, in order the keep the extent of the wave- 
function much smaller than the extent of the field, a too must be kept large 
enough. 



6 Connection with the adiabatic approxima- 
tion. 



It is instructive to study the same problem in the limit where the system is 
extremely adiabatic. For simplicity wc specify to the case G = 0. In this 
approximation, the direction of the spin n is locked to the direction of the 
local magnetic field, so that 



n 



H Jl^z 



H Vlyi 



(25) 
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Substitution of Eq. 
Eq.(RD gives 



into Eq.(|^ and discarding the equation for the spin 



m 



— + 2fi^z X — + n;zx 
dV at 







S 


X r) 


= -/iV 


H VLj-z 




/i 



(26) 



with H given by Eq.(|). 

It is important to note that this approximation is different from the 
time-averaged orbiting potential (TOP) approximation. In the latter, one 
constructs a time-dependent potential Vxapij^'t) oc |H(r,t)|, then averages 
Vtop(i',^) over time to get a time-independent potential V^Qp (r). Here, we 
work in the rotating frame, in which H (r) is time-independent, and construct 
a time-independent potential 



VAD{r) oc 



S 

H 

/i 



Normalizing the Cartesian components of the position vector r = xx + 
yy + zz to Rq, we find that a stationary solution to Eq.(^) is given by 



1 



2fi2 
zq = Q. 



(27) 



This result agrees with Eqs.(|Tl]) and ([T2|) for the case g = 0. The other 
possible solution, in which xq is positive, is discarded because it is not stable. 
Substituting 



X = Xq + 5x 
y = yo + Sy 

Z = Zq + 6z 



(28) 



into the normalized form of Eq. (^6]) , and expanding to first order in the 
perturbations gives 



d^x _ ^^dSy _ ^ Q 



df^ " dt 



dt' 



dt 
d'5z 



4fi2a + 1 
-6z. 



5y 



(29) 



a 
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We therefore find that the axial translational degree of freedom is decoupled 
from the rest, with a frequency 



which agrees with Eq. (|2^) . For the lateral translational degrees of freedom 
we have 



6x\ / 



-2iiun -uj^ - + — \ J V 

AQ^a + 1 



M 

for which a non-trivial solution exists whenever 

{An^a + 1) det M = {Ail^a + l) lu^ - (Sfi^a + 3^^) uj"^ + 4^]6^ = o. 
This equation furnishes the frequencies 



2 (4fi2a + 1) 

corresponding to the laboratory-frame frequencies 

V oa 

This result also agrees with Eq . (|23|) . 

Though in this model the spin is locked to the direction of the field, we 
can nevertheless define a precessional frequency Up by calculating the field 
at the stationary point Hg and define ujpQo = l^Hg/S. This gives 

1 

ujp = a + 



4fi2 



which for aVL^ ^ 1 coincides with Eq. 

Note however, that the adiabatic approximation presented here holds 
whenever the precession speed ujp is large compared to the vibrational fre- 
quencies. In addition, ujp should also be large compared to the rotation 
frequency Vt. The extreme case ujp = VL defines a line in the a-Vt plane which 
for large a approaches asymptotically to the line a = VL. 
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7 Discussion 



We have shown that our exact analytic results reduce to the formulae derived 
via the time-averaged orbiting potential approximation, in the case where a 
is large. In addition, the stability diagram that we found, suggests that the 
TOP trap is very flexible for the experimentalist in terms of the allowed 
parameters. We have also shown that under the adiabatic approximation, 
where the direction of the spin n is locked to the direction of the field in 
the rotating frame, we recover, for large magnetic field a, the exact mode 
frequencies. 

It is interesting to note that Eqs. (|^ and pave the way for a quantum- 
mechanical treatment of the same problem in the comoving frame. When 
gravity is neglected, the Hamiltonian for this system is given by 

2m \ fi J 2 ^ ' 

where A is a vector potential field satisfying 

V X A OC i^r Z, 



is the spin S (where S could be either or 1/2 or 1 etc.) vector of Pauli 

S 

matrices, P is the vector momentum operator, H fi^ is the magnetic field 

as seen in the rotating frame, and — mf2^ (x^ + y^) /2 is the centrifugal po- 
tential. In order to diagonalize the magnetic part of the Hamiltonian, one 
performs a local passive transformation of coordinates on the wave function, 

such that the spinor is expressed in a new coordinate system whose 2;-axis co- 

S 

incides with the direction of the local magnetic field H at the point r. 

/i 

This rotation does not affect either the centrifugal term or A. The momen- 
tum P however, transforms to P — A' {v,a^, ay,a^), where A' (r,(Tf , cry,a^) 
contains non- diagonal elements as it includes the spin degree of freedom. For 
typical values of parameters, the non-diagonal part may be treated as a small 
perturbation, and the lifetime of the particle in the trap may be calculated. 
This technique has already been applied successfully to a ID toy-model time- 
independent magnetic trap[^, and to a loffe-like 2D trapfl^, for the case of 
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spin S — 1/2 particles. As the TOP trap is used to capture Bosons, it is 
more resonable to study it for the case S = 1. We beheve that despite the 
complexity that the Coriolis and centrifugal forces add to the problem, it 
is possible to solve TOP trap quantum-mechanically. The analysis of this 
problem is still under study. 
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Table 1: Comparison of mode frequencies 



Mode 


Measured 


TOP [| 


Exact Anlaysis 


Prec. freq. MHz 
Axial freq. Hz 

Lateral freq. Hz 


67 ± 1 
24 ± 1 


~ 7 
69 ±2 

24 ± 1 


7 

6 

] 


.38 
7.99 

f ±24.44 
15ir± 24.44 



Figure 1: Stable region for (7 = in the a-Vt plane. 
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